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I. INTRODUCTION 



The standard definitions of intensive tliermodynamic parameters, such as temperature, 
seem to require the system in question to be in thermodynamic equilibrium. In this paper, 
we explore the entropy and temperature of radiation out of equihbrium and show that, 
within hmited restrictions which do not require equihbrium, the radiation temperature is 
well-defined and distinct from any associated matter temperature and from plausible but 
incorrectly applied equilibrium definitions. We use laser radiation as our major example. 

Laser radiation is a fascinating example of a highly organized quantum system of quasi- 
coherent bosons. ^'^'^ A laser beam is supported by external pumping, which keeps the beam 
far from thermodynamic equilibrium. The laser shares this feature with other steady-state 
systems that are kept from equilibrating by external constraints. We find that the tempera- 
ture of laser radiation far exceeds the temperatures of the laser cavity and the lasing atomic 
transition. Misidentifying the radiation with the matter temperature leads to erroneous 
estimates of the laser radiation temperature that are as much as ten orders of magnitude 
too small. 

Photon number, unlike energy, is not conserved. The Gibbs-Duhem relation for radia- 
tion, SdT — VdP = 0, implies that the two intensive thermodynamic parameters, pressure 
P (conjugate to volume) and temperature T (conjugate to energy), reduce to one inde- 
pendent intensive parameter, which is usually identified as T. This feature of radiation 
thermodynamics, like the photon's zero mass and lack of rest frame, makes radiation ther- 
modynamics much simpler than that of matter, which has conserved particle numbers and 
nonzero chemical potentials.i^ It also makes generalizing intensive thermodynamic parame- 
ters out of equilibrium much easier. Thus radiation is a natural context in which to introduce 
nonequilibrium temperature. 

A properly defined nonequilibrium temperature has physical meaning. It occurs in the 
entropy production rate S, an important measure of both how far a system is from equi- 
librium and how fast it is approaching equilibrium.^ S has a generic form rooted in the 
equilibrium expression for the entropy differential, dS = dQ/T, where dQ is the differential 
of heat or random energy, the change in system energy while holding volume and particle 
number constant. The form S ~ Jq(12)(1/Ti — I/T2) expresses the entropy produced by 
two subsystems (1,2) at temperatures Ti, T2 as they exchange a heat flux Jq(12). Subsystem 
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temperatures occur naturally in expressions for entropy production. S is positive semidef- 
inite and vanishes if and only if Ti = T2, the condition for thermal equilibrium. The heat 
flux Jq(12) vanishes in this well. 

The difference 1/Ti — I/T2 of the inverse temperatures is a measure of how far out of 
equilibrium the two subsystems are. The heat flux Jq(12) is a measure of how fast the sub- 
systems are approaching equilibrium with one another, assuming no external pumping of the 
system. (With external pumping, Jq(12) is a measure of how much power has to be injected 
into the system to keep it from equilibrating.) The product of these two quantities, given 
by S, combines the two measures into a single quantity characteristic of a nonequilibrium 
process. 

The books by Reichl^ and De Groot and Mazur- explain in detail the significance and 
role of entropy production in nonequilibrium matter systems. Section IVIII and Refs. and 
explore entropy production in radiation and radiation-matter systems. 



II. MISAPPLICATIONS OF TEMPERATURE OUT OF EQUILIBRIUM 

Many formulas for various physical quantities have units of temperature. These formulas 
frequently are defined as the temperatures that equilibrium systems would have if the energy, 
entropy, number, etc., were all rearranged in some particular way. There are many ways 
to rearrange a system into an equilibrium state. Consequently, many definitions of these 
pseudo temperatures are possible, but none represent a temperature of the actual system, 
and the application of temperature to nonequilibrium systems is ambiguous.— 

Consider a helium-neon laser with the transition line Aq = 632.8 nm. A heuristic but 
incorrect equipartition argument which sets the mean photon energy equal to ksT would 
associate this transition with a temperature ksTo = Hc/Xq or Tq ~ 2 x 10^ K. The quantity 
To is the temperature that a black body would have if the distribution of energy among the 
photons were contrived to have an average photon energy corresponding to the photons in 
question. But the rest of the distribution is far from a blackbody for any laser. Tq refers to 
a photon gas with a different distribution, energy, and entropy. So Tq has no direct physical 
meaning for laser radiation, which is far from equilibrium with itself and with its matter 
source. 

Laser radiation often is idealized as having infinite temperature, as lasers are interpreted 
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as a source of pure work, although ultimately this interpretation is unsatisfactory. It would 
put a powerful X-ray laser on an equal footing with a pocket red diode laser pointer, powered 
by watch batteries. Is there no difference in "temperature" ? There seems to be none if both 
lasers are described by infinite temperature. 

Other possible pseudo temperatures include the temperatures that a beam would have if 
the same energy or the same entropy were arranged differently, for example, in a black body 
distribution. Such temperatures do not reflect the actual distribution of energy, entropy, 
and photon number in the beam. Different as these definitions are from each other, they 
would all agree if the laser radiation were forced into thermodynamic equilibrium. 

There are other nonequilibrium temperature definitions that reflect the actual energy and 
entropy distributions. For example, objects in the laboratory surrounding the laser have 
meaningful local temperatures. The material around a He-Ne laser is at room temperature 
Troom ^ 300 K, and the He-Ne gas is at T^^ ~ 400 K. 

The inverted populations N2 and A^^i of the upper and lower laser energy levels E2 and 
El are associated with a temperature T21 through a formal Boltzmann distribution, 

N2/N1 = exp[-{E2 - Ei)/kBT2i], (1) 

which does not hold for all levels. Such a pseudo temperature can be defined for any two 
levels. In this case, the definition and population inversion imply a negative value for T21. 
But such a definition lacks a thermodynamic justification. 

Although an inverted atomic population is essential to lasing action, there is no reason 
to attribute the temperature T21 to the radiation field, a separate entity with its own ther- 
modynamics. In Sees. IIVI and |3 we will identify the natural radiation temperature that 
represents its distinctive nature, and that points to a more general definition in terms of an 
energy distributed over states. Temperature emerges not as a proxy for energy, but instead 
as a measure of how the energy is organized among the various microstates. It possesses 
a rigorous definition and plays a natural role in nonequilibrium systems that go beyond ad 
hoc heuristic estimates.— 
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III. DYNAMICS OF THE RADIATION FIELD 



We analyze the radiation field in a box and decompose it into plane wave modes.i^ 
Each field mode of wavevector k and given polarization fills space. (We ignore polarization 
in this paper for simplicity, which does not change the generality of the argument.) The 
fundamental mode variable is its complex amplitude = OkC*'^'', in terms of its modulus 
Ok and phase 0k- 

Each mode is a linear harmonic oscillator and can be mapped onto a quantum harmonic 
oscillator in terms of lowering (raising) operators Sk (Sj^), satisfying [ak,aj^] = Ik, for each 
mode. The state of a mode lives in an infinite-dimensional Hilbert space. The occupation 
number basis \n) forms a complete orthonormal set of number operator eigenvectors corre- 
sponding to different photon numbers: ri-kl^k) = fikfikl'^k) = '^kl^k), and X]nl'^)(^l ~ ^■ 
The field Hamiltonian H = "^y. hv{hy_ + 1/2), with v = ck. The state of the whole field is 
the direct product of all its mode states and lives in a field Hilbert space (Fock space). 

An alternative basis is the set of coherent states |ak), eigenstates of the lowering operator: 
a^kl^k) = ttkl^k); defined so that the occupation number expectation A^k = (o^klflkOklctk) = 
l^kp = Ok- Coherent states satisfy 



The state |ak) is the quantum analogue of a single classical mode k with amplitude a^. 
Physical results can be obtained with either the coherent state basis or the photon occupation 
number basis |n).— 

Problem 1: Use the definition of {n\a) to derive the identities above: the result for {(3\a) 
and the integral identity for the coherent states \a). Also show directly, by expanding \a) 
in In), that d\a) = a\a). 

Thermodynamics requires a statistical description with an ensemble of many Fock spaces. 
Each mode has an ensemble of amplitudes and phases, and the ensemble of the whole field 
is the direct product of the mode ensembles.— 

Very high temperatures are to be expected in low entropy, high energy bosonic systems 
far from thermodynamic equilibrium such as the laser. Such high temperatures are natural 
for bosons, because they lack an exclusion principle. Ideally, a multiparticle bosonic system 
could be driven to zero entropy by putting all particles in one single-particle state, without 





5 



changing the system's total energy. As shown in Sec. IIV CI if the total energy remains fixed, 
the temperature of that state would be infinite. 

IV. NONEQUILIBRIUM THERMODYNAMICS 

A system's statistical ensemble defines the system's entropy. Defining the temperature 
of a system or subsystem necessitates the restriction to cases where its entropy and energy 
have a functional relationship. A physically meaningful definition of temperature should 
refiect the system's actual state and properly reduce to equilibrium temperature. 

A. The Entropy in General 

To define the entropy, we introduce an ensemble of a large number M of copies of a 
system assumed to be made up of discrete, countable "things." Physically, these things are 
the system's fundamental degrees of freedom. Because of statistical fiuctuations, each copy 
is microscopically different and distinguishable from the others. 

Consider the entropy of two system copies. Let Wi be the overall number of ways that 
the first group of things can be arranged, and W2 the corresponding number for the second 
group of things. The overall number of ways W that the two groups can be arranged is the 
product W = Wi X W2- For the entropy, we want a function of W that matches the usual 
definition for the special case of equilibrium. Therefore, this function of W must depend on 
the two combined subsystems by adding the entropies of the separate subsystems. Such a 
function must be a linear function of the logarithm of W. The definition S = /cslnVT, up 
to a multiplicative factor and an irrelevant additive constant, is the only function that does 
so. Note that entropy should be additive, but need not scale in a simple way with the size 
of the system. For a nonequilibrium system, although subsystems contribute additively to 
the whole, the subsystems do not have to be homogeneous and contribute to the whole in a 
simple scaling fashion. 

Each system copy has the same internal probabilities for being in any particular state 
a. The number of system copies in state a is Mp„ = m^. Suppose that the ensemble has 
mi systems in state ai, m2 systems in state (T2, etc. The number of ways this ensemble can 
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be created from the system copies is 



WM = M\/[mi\m2\...mJ...l (3) 

and the entropy for each member of the ensemble is S* = ks^'^iyVM) / M . For M large, 
Stirling's approximation yields 

MS = kB[M\nM -^{m^lnm^)]. (4) 

(7 

(Stirling's approximation is n\ ~ vJ"- \J l-nn for n large. We ignore additive constants in 
the entropy, as these are independent of the system's thermodynamic state.) 
Because rrif, = Mp^-, 

S = /cspnM- J](p,lnM + p,lnp,)] 

(7 

= fcB[(l-5^Pa)lnM-^p,lnp,] 

(7 a 

= -ke^Pa'^^Pa, (5) 

cr 

where the last line follows because "^^Pa = 1 and the ensemble's system copies are inde- 
pendent. 

A thermodynamic system's statistical ensemble is equivalent to its normalized, Hermitian 
density operator p : Tt{p) = Xlo-Po- = 1- The p^ are the eigenvalues of p, with eigenvectors 
I a). Any observable O has a statistical average 

{0) = TT{pd) = Y,PMd\a). (6) 

(T 

The average of the operator — fc^lnp is the entropy, S = — A;BTr(/31np). The mean energy 
E = Tt{pH) is the average of the Hamiltonian if.— 

B. The Validity of Nonequilibrium Temperature 

Out of equilibrium, the entropy 5* lacks a clear functional dependence on the total energy 
E, and the definition of T becomes ambiguous. However, an important generalization of T 
is still possible, with certain restrictions. The crucial requirement is that the whole system 
decompose into subsystems whose entropies are functions of each subsystem's energy only. 



7 



That is, the nonequihbrium system must decompose into subsystems each with its own 
equihbrium. Restriction to full equilibrium is unnecessary. 

We write the complete system Hamiltonian as H = J2p^P + J2q>p^pg^ where the first 
sum is over separate subsystems and the second is over all interactions between subsystems. 
The full p evolves via the quantum Liouville equation ihdp/dt = [H-p]- For a particular 
subsystem p to have well-defined thermodynamics, the necessary restrictions are the follow- 
ing: 

(1) The whole system's density operator p must be factorizable into independent subsys- 
tems: p — ®pPp. This factorization implies that the total system entropy S is a sum over 
subsystems: S = Sp. 

(2) The commutator [Hp, pp] must be negligible or zero. Then pp and Hp are simultane- 
ously diagonalizable, and a functional relationship pp — Pp{Hp) is possible. The subsystem's 
ensemble probabilities are — {a\pp{Hp)\a) — pp^{Ep^). That is, each eigenvalue Pp^ of 
Pp is a function of the corresponding eigenvalue Ep^ only. It is not a function of the other 
energy eigenvalues Ep^ of subsystem p or of the {Eg^} of the other subsystems q. (If the 
interactions Hpq between subsystems are negligible, then the total system energy £" is a 
sum over the Ep^s and each subsystem is in a stationary state. But this restriction is not 
necessary. For example, laser radiation modes in the laser cavity are strongly coupled to 
external pumping.) 

These special conditions allow a functional relationship between ensemble averages of 
entropy and energy for each subsystem that does not depend on the specific probability 
distribution {pa}- That is, Sp = Sp{Ep), which leads to a natural temperature of subsystem 
p: 

1 = (7) 
Tp dEp 

where 7], is a rate of change between extensive quantities, as it normally is in thermodynam- 
ics. Such a relationship exists for photons if the subsystems are chosen over small enough 
ranges of energy. We note that (i) the functional relationship Sp = Sp{Ep) is an equation of 
state, (ii) If Sp is a function of Eg, for q ^ p, the temperature Tp is still defined, but is not 
a function of subsystem p alone, (iii) If many subsystems are in contact with one another, 
their temperatures represent a basis for their complete thermal equilibrium. For each pair 
of subsystems (p, q), Tp — Tg — T, where T is the single temperature of the whole system: 
l/T^dS/dE. 



A familiar example of such a subsystem temperature occurs in the local thermodynamic 
equilibrium of the radiation field typical of stellar interiors. In this case, periodic boundary 
conditions are applied in small local boxes at positions r. A different field for each Fourier 
mode k, ensemble, and temperature T(k) can be defined in each box. (Very long wavelength 
modes not fitting into the small local boxes must be ignored.) Thus the subsystem is a 
particular mode k in a box at r, and Tp becomes a photon phase space temperature T(r, k). 
At a point r in space, this brightness temperature of the photons is a function of mode 
frequency u = ck and direction k and is conventionally written as T^{^). In Sees. IVland lVH 
we show how this temperature emerges naturally as a genuine thermodynamic one. 

C. Pure States and Infinite Temperature 

A pure state has zero entropy. One of the p„ (say p^) is unity, while the others vanish, 
and so 5* = 0. There are many ways to realize a pure state with an idealized laser. The two 
simplest are a pure coherent state p = or pure occupation number state p = \n){n\. 

Problem 2: Generally, p is a mixture of projection operators over different states; for 
example, in the number basis, p = '^Pnm\n){m\. But the ensemble of a pure state is 
p = a single projection operator, not a mixture. Prove that the eigenvalues of a 

pure ensemble can only be one or zero. Hint: A normalized projection operator 11 = 
with = 1, satisfies 11^ = II. Then infer that the eigenvalues p of a pure p satisfy the 

formally identical equation, = p. 

We define a single mode's ensemble {Pa}- The ensemble average energy is given by 
E = J2aP(^^^' where E^j = {a\H\o-). The subsystem temperature is defined by a derivative. 



where dp^j = 0. 

This definition of subsystem temperature is a simple restatement of the thermal equilib- 
rium of subsystems across ensemble members. This argument does not require all members 
of the ensemble to participate, but the probabilities are assumed to be normalized over the 
active members of the ensemble, with the rest are ignored. Equation (jHl) thus provides a 
legitimate equilibrated temperature of only those ensemble members indicated by the sum- 
mation. 
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Y,^E„dp„ ' 



(8) 
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Consider two cases based on a sequence of equilibrated subsystem probability distribu- 
tions with the temperature defined by Eq. (jH)): 

(1) In equilibrium across all ensemble members, the p„ are very small, as the probability 
is spread over the entire ensemble. As the mean ensemble energy E falls, while maintaining 
equilibrium, the average is dominated by the lowest-lying value of the set {-E'er} as the lower 
and lower energy states are occupied. Then {pa} becomes dominated by one pr ^ 1, while 
the others vanish. But in most systems, the denominator in Eq. (jHJ vanishes faster than the 
slowly-changing logarithm in the numerator. So in equilibrium, T — ^ as -E and S ^ 0. 
The classic example is the Bose-Einstein condensate.— 

(2) In the nonequilibrium case of interest in this paper, the mean ensemble energy E is 
fixed. Now redistribute this energy and the ensemble probability into ensembles restricted 
to successively smaller subspaces of the Hilbert space. Equilibrate the energy on these 
successively smaller ensembles. Equation (jSI) still holds on these smaller ensembles, but 
E and the {p„} are distributed over fewer states. One p^ approaches unity and the others 
vanish. The number of ensemble members declines toward one, the pure state limit, yielding, 
in the limit of one possible subsystem in the ensemble, the limit 

T^~-^-0 (9) 

for the temperature. So T oo for a pure state if the total energy is fixed while the entropy 
vanishes. 

These results show how temperature is related not only to the ensemble energies E„, but 
to the ensemble distribution {pa} as well. 

The radiation field of a real laser is not in a pure state. Each mode has a statistical 
ensemble of amplitudes or occupation numbers, with an associated entropy. Thus, whether 
the field is given a classical or quantum description, a real laser does not have infinite 
temperature, although laser temperatures are often idealized as infinite. 



V. RADIATION ENTROPY AND TEMPERATURE 

From this point, we will work with individual photon modes k and drop mode labels where 
they are not needed. Statistical ensembles of radiation can be described by density operators 
p.i^ Although a real laser is not in a coherent state, the coherent state basis is exceptionally 
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useful, because it connects the radiation mode description in terms of amplitude and phase 
with photon number and energy. 

In principle, we could calculate the entropy S and other thermodynamic functions from p 
and find T. But S is difficult to calculate for arbitrary p. We consider only the special case 
where the radiation field's statistical ensemble is independent of the phase 0. For most lasers, 
the phase is fully randomized. The mode energy E = hvN is always phase- independent, 
depending only on the ensemble's photon number expectation = Tr (a^a). Thus the phase 
does not affect or T in this restricted case. This restriction simplifies S{N) to a general 
nonequilibrium form which we can obtain by a simple counting argument appropriate for 
bosons.— 

A. Coherent State Basis — Random Field Phase 

The density operator can be expanded over coherent states, p = J(faP{a)\a){a\, and 
then projected on to the \n) basis, Pnm = {n\p\m). The normalization is fixed, Tr(/)) = 
J£a P(a) = ^^p„„ = l.i^ 

The entropy is given by S* = — A;BTr(plnp) = — A;^ ^^(n|pln/)|n). We restrict the dis- 
tribution function P{a) to be independent of phase (p and depend only on the modulus a. 
The matrix pnm then becomes diagonal: 

//■oo ^271+1 ^ 
(fa{n\a){a\m)P{a) = 27i6nm da — — e"" P(a). (10) 
Jo 

The entropy simplifies to 5* = — fc^^^pnlnp^. The eigenstates \a) of p are the occupation 
number states |n). The distribution P{a) is phase-independent if the phase distribution for 
the amplitude a is random over the interval G [0, 27r). 

The restriction to random phase and the result that the mode state can be character- 
ized by counting the photon number alone validates the phase space approach to photon 
thermodynamics. For each mode of the field labeled by a wavenumber k, the mode state is 
determined solely by its mean photon occupation number, = '^nPk„ny_M 

B. Phase Space Reduction: Counting Photons 

Making the field's statistical ensemble phase- independent allows us to find S{N) by simply 
counting states for identical bosons, treating the field quanta as particles. In this case, all 
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of the entropy is due to the randomness of the phase. The reasoning of this section parallels 
that of Sec. WJ\ 

Consider N identical things and G identical possible places to put them. Imagine not 
knowing where the things are among the G places. Then there are (A^ + G — 1)!/A^!(G — 
1)\ = W ways of arranging the things among the G places. Now introduce an ensemble 
of M system copies, with each copy labeled by a. For systems such as lasers with large 
numbers of photons, we can assume that G and ^ 1 and use Stirling's approximation. 
G and A^ depend on the system copy a in the ensemble. The entropy is then: 

MS = A;B^[-G,lnG^ + (G^ + Ar,)ln(G^ + Ar,)-Ar,lnAr,] 

cr 

= kB^G4{l+n)ln{l + n) -nlnn], (11) 

where n„ = N„/G^ is the mean occupation number in the ensemble and approaches the 
cr-independent limit n = N/G, with G = as M — >• oo. Because hu/T = dS/dN 

and A^ = Gn, 

k^T = , , , (12) 
ln(l + l/n) ^ ^ 

Problem 3: Combine Eqs. (Illj) and ()12j) to express S in terms of T. 

To connect the result ()12|) to radiation observables, identify the G "places" with 
the Fourier space of modes k and use the general relation of mean occupation number 
n to the specific intensity '■ n = c^/;/(k)//ii/^, where is the radiation en- 

ergy/area/time/frequency/solid angle.-^^^ The resulting expression for the temperature, 

kBTM = — — ^, (13) 

can also be derived from the Planck blackbody expressions for ly by solving for T and 
defining the resulting temperature in terms of ly. The value of ly can then be arbitrary.— 
The radiation temperature varies as a function of beam direction k as well as frequency u. 
Except at zero frequency, Ty is zero only for vanishing n or ly (photon vacuum). 



C. Temperature as a Lagrange Multiplier 

Instead of Eq. (jH)) with its assumptions about probabilities, the temperature can be con- 
sidered as a Lagrange multiplier for holding the ensemble mean energy fixed while maximiz- 
ing the entropy. This procedure can be extended to subsystems, such as a single radiation 
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mode or a collection of modes constituting a Fourier subspace K, a part of the whole mode 
space. The conclusions of the previous sections are confirmed in a different way, as long as 
the mode phases are random. 

If we express Eq. (II 1|) in Fourier space, the real space volume density s of entropy con- 
tributed by photons of wavenumber k and frequency u = ck in the Fourier subspace K 
becomes: 

f d^k 

s= kB[{l+n)\n{l + n)-n\nn]-——, (14) 

where n is the mean occupation number. The real space volume density e of energy in K is 

d^k 
(2^3 



nhv . (15) 



'K 

The maximum entropy at fixed energy is given by 



= 6(8 + Pe) = (i±^) - Phv]-^^, (16) 

from which we conclude that 

We find an apparent black body distribution for the specific intensity: 

hv^ 1 



(18) 



and identify the Lagrange multiplier (3 = l/ksT. We define the specific intensity of the 
entropy radiation J,^ by photon counting, the same way that the specific energy intensity 
is defined.— Then we find another version of Eq. (|7|): 

^ = ^ (19) 

T,(k) diM 

If K were all of Fourier space, this distribution would be an equilibrium one with a 
single temperature T. But we have not restricted K. If K is any part of Fourier space, 
the distribution is a black body distribution, but truncated to only its domain of Fourier 
subspace. This subspace is complete with its own temperature kBT[K] = 6e[K]/6s[K], even 
though the whole distribution is not present. The simplicity of radiation thermodynamics 
converts the functions s = s{n[K]) and e = e(?T,[i^]) into an implicit function s{e) over 
K, matching the more general argument of Sec. IIVBI Therefore each K forms a distinct 
thermodynamic system with a temperature given by Eq. (fT3j) or obtainable by inverting 
Eq. (UHl). 
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This argument can be extended to a of zero volume. By squeezing a fixed radiation 
energy into a vanishingly small frequency interval Az/ and beam solid angle cone AQ, we 
recover the infinite temperature of a pure state for one mode k, with u = ck, as we now 
show. 

Consider a fixed arbitrary specific intensity distribution and break it into two parts, 
a black body function B^, of constant temperature T, plus a finite deviation distribution 
Dy(k), which may depend on the direction k. Thus = + Du. The energy density u 
within (Az/, Af2) is given by 



cu 



/ dn^ / dul^ = dn^ dv{B^ + 
Jaq Jau J An J Au 



dfik / dv ^l^lfjf + [ dn^f dv D,. (20) 

If we change to the variable x = hu/ksT, the first integral of the right-hand side of Eq. (j20| 
can be recast as: 



rfa/ dx-^. (21) 

It is clear that letting Au go to zero implies that Ax goes to zero also. From Eq. T 
cannot reach zero except for zero intensity, — ^ 0. This limit cannot be reached because we 
require that this beam have fixed, finite energy density u. Because the integrand is finite, 
the double integral in Eq. (j?T|) vanishes in the limit of Au and AQ going to zero. 

It follows that the second double integral on the right side of Eq. (I20j) vanishes in this 
limit, as does the double integral in Eq. (j2T| . If we use Eq. (j2T| to rewrite Eq. (j20|l . we 
obtain 



lim 



, , \im I (cu-[ dn^ ! duD^)\ [ dVl^ [ dx^^] ^1 = cx). (22) 

Thus T becomes infinite when a fixed energy is concentrated into a single mode. This limit 
recovers the special ensemble of a pure state discussed in Sec. IIV CI fixed energy with zero 
entropy. The ensemble degenerates into a single Hilbert space, with an exact photon number 
A^k = E/hu or field mode modulus Ok = \J E/hu. In this limit, these ensemble averages 
have no statistical uncertainty. The mode fills all of real space and has one wavelength and 
direction, a fixed absolute amplitude and energy density, but random phase. The phase 
is independent of energy and does not affect T. This result formally justifies the idealized 
picture of laser radiation as having infinite temperature. 
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When radiation is in equilibrium with matter, all Fourier subspace temperatures become 
the same finite value, returning Ii, to a full black body distribution, as in Eq. (fTHj) . but 
valid over all modes. If the matter is selective in its frequency or directional response to 
radiation, the relaxation of the radiation to equilibrium is similarly limited. (The radia- 
tion intensity distribution can be found in general by solving the equation of transfer that 
describes radiation transport through matter.—) 

D. Temperature versus Intensity — Classical Limits 

For a fixed and frequency-independent ly = I, very different temperatures are found at 
different frequencies. We rewrite Eq. fjl2|) and ()13|) in terms of the reciprocal of the mean 
occupation number, z = hv^ /{c^I^) = l/ra, 

kBn = h(^^y'g{z), (23) 

where 

g{z) = z'^/\n{l + z). (24) 

The gain function g{z) is plotted in Fig. ^and determines how the intensity and frequency 
are related to the temperature for radiation. g{z) is singular at 2; = 0, with a minimum at 
approximately z = 15.8; it grows gradually, unbounded, with increasing z. 

Clearly "hotter" radiation sources for a given energy are either at low or very high fre- 
quencies. For example, a laser in the 600 nm range is less "hot," watt for watt, than an 
X-ray laser with wavelengths of the order of Angstroms. This result makes sense because 
Eq. (|13|) must give a constant value for a black body, having a minimum in the middle 
of the frequency domain. Alternatively, the curve represents points where the black body 
distribution crosses a constant specific intensity /. It crosses at two frequencies except when 
the temperature falls so low that it does not cross at all. 

The low- and high-frequency limits of Eq. (PHj) can be understood differently. If we count 
photons as particles in terms of their energies e = hu, the relation n = [exp(e//cBT) — 1]~^ 
is independent of h. In the limit e^ooorT-^0, n— > exp(— e/fc^T), the Maxwell- 
Boltzmann distribution. Photons in this limit act as classical particles, and n tends to be 
small {z — > 00). Then consider the low-frequency or high-intensity limit of Eq. (j^Hj) . In 
terms of intensity, the resulting classical relation ksTi, = (c^/z/^)/ is independent of h and 
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arises from a set of classical thermal oscillators. In this Rayleigh- Jeans limit, z ^ and 
n ^ oo. 

Photons therefore have two different classical limits, at high frequencies as particles and 
at low frequencies as a classical field. The field limit is possible because photons are bosons, 
and large numbers of photons can coexist in the same field mode. Coherent states can be 
constructed as analogues to classical field states. 

VI. THE TEMPERATURE OF A REAL LASER 

Although lasers do not have infinite temperatures, even common, low-power lasers have 
temperatures closer to those of stellar interiors than to everyday matter temperatures, far 
exceeding the mistaken equipartition estimates of Sec. m Lasers operate far into the high- 
intensity Rayleigh- Jeans limit of Eqs. (fT^ - (fTn|) . where ksT^, = {c^/u'^)!^. To see these very 
high temperatures, the specific intensity needs to be extracted from the laser power P. 

The flux density F is related to the specific intensity Iu(^) by 



the wavevector and surface area normal vector f . is the mode solid angle element, and 
dA represents the differential exit aperture area of the laser beam. 



from each point of the aperture would be constant over the aperture, F = P/A. The power 
would be 




(25) 



The surface area element dA = r 



"^dQr, and cos^ = f ■ k is the cosine of the angle between 



If /i/(k) were constant with beam angle and frequency, then /i/(k) = J, and the power 



P = vrA/sin^ ^1/2, 



(26) 



where 61/2 is the beam spread half-angle, and A is the aperture area of the laser. 
Let the specific intensity have the factorized form: 



/,(k) = Fo$(i/)P(k). 



(27) 



Each function, $(i^) and T>{'k), is separately normalized to unity: 




(28) 
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The solid angle variables are 9 and (p. They are the polar angle relative to, and the azimuthal 
angle about, the beam axis, respectively. 

We will consider typical lineshapes and angular distributions.— The simplest lineshape is 
the Gaussian form arising from Doppler broadening by thermal agitation of the lasing gas, 

= exp[-(z/ - vof/{AvDf]/^/\vD, (29) 

and Az/£) is the half-width due to Doppler broadening: 



^i^D = ^oV (30) 
V mc^ 

m and Tgas are the mass and temperature of the gas atoms emitting the radiation. 

The simplest angular distribution is also a Gaussian from the laser cavity resonating in 
its fundamental mode. 

V{e) = 2 exp[-2eV^o]/^^0' (31a) 
j dipj d{cos9)coseV{9) = 1. (31b) 

The distribution is azimuthally symmetric about the forward beam direction, with 6o as 
the half-angle, e~^-power, beam divergence. An ideal laser's beam divergence is diffraction- 
limited at the aperture: 6q = 2Xq/ttD, where D equal to the aperture diameter. 

Our results are for a red He-Ne gas laser, wavelength Aq = 632.8 nm or uq = 4.741 x 
10^^ Hz. We assume a power of Pq = ImW. Figure |21 is a plot of the radiation temperature 
T{v^9) for D = 1mm and Au^, = 0.9 GHz, corresponding to Tgas — 390 K. The beam 
aperture area A is ttD'^/A ~ 8 x 10~^m^, and the total flux density Fq = Pq/A ~ 1.3 x 
lO^W/m^. The beam divergence Oq ~ 0.4mrad. 

Figure El shows the radiation temperature T(i/ = z/q, 6) for several beam diameters D and 
corresponding beam divergences 6q, with the same Pq and Aud. The peak temperature is 

= UQ,9 = 0)= . (32) 

It is independent of D, because the normalization of temperature is independent of D: 

T ~ Fo/el ~ Po^o"V^ ~ e^^D-^ ~ D^D-\ (33) 

Only the shape of the angular distribution depends on D. Note that T^, diverges if A^) or 
6*0 0, reproducing the pure state of Sec. IV Gl 
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To infer the radiation temperature of a laser requires knowing the specific intensity I^, of 
its hght. Because a laser's radiation is so intense, we can use the Rayleigh- Jeans limit of 
Eq. (fT^: ksTi, = (c^/u'^)!^. Finding I^, requires making separate geometry- and frequency- 
based measurements. The former requires an intensity measurement at a fixed frequency, 
typically with a photodetector. The latter determines the spectral characteristics of the 



We place a photodetector of known area A and frequency response or measurement 
efficiency -R(z^) orthogonally across the beam and centered along the beam axis {6 = 0), 
at a distance r from the beam exit. -R(z^) is normalized so that its integral is unity. The 
photodetector measures an incoming power P: 



with the definition ()25|1 . The factor of 2tt represents azimuthal integration around the beam 
axis and assumes axial symmetry of the beam's intensity. 
We make the following simplifying assumptions. 

1. Assume that the frequency and angular profiles are Gaussian, like the expressions 
()29|) and (jHT|) . From these profiles, we can infer the peak fiux density Fq with a few 
measurements. 

2. The solid angle subtended by the photodetector is f2 = A/Anr^. Assume the photode- 
tector itself is circular. The half-angle 6a it subtends with the beam axis is given by 
tan 9 a = ^/A/7^r^. 

3. All relevant angles, 6a and ^o, are small. Then cos 6' ~ 1 — 6'^/2. The angular integra- 
tion becomes: 



4. The frequency integration can be simplified if we assume that -R(z^) is constant across 
a range AF, R{i') = (AF)^^. If we do the frequency integration numerically, however, 
this simplification is not necessary. 

We can use values of z/q, Aud, and 6q that are specified by the laser's manufacturer. With 
additional instruments, we can measure these quantities, although such measurements are 



beam.' 



2.16.17 




(34) 




(35) 
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more difficult than the power measurement. For example, we can determine the frequency 
profile and infer uq and Auo by measuring the specific intensity as a function of wavelength. 
Light wavelength is typically measured with a grating spectrometer. But a grating does 
not have the frequency resolution needed for a laser, and instead we should use a Fabry- 
Perot interferometer-based analyzer. We can determine the angular profile and infer 6q by 
measuring how the relative intensity varies as we vary r, the distance of the photodetector 
from the beam aperture. Because we move the photodetector farther away from the laser, 
it subtends a larger half-angle 6a with the beam axis. 

If we substitute these derived values and assumed simplifications into Eq. (jH^ . we can 
do the angular integral analytically. The frequency integral can be done numerically. With 
the measured P, A, z/q, Aud, and 9o, we can infer Fq and /^(z/, 0) from Eqs. (j27j), (j^ . (jSH), 
and ()34|) . Your results should be similar to Figs. |21 and El unless your laser is not operating 
in its fundamental frequency mode. (If the photodetector is calibrated in absolute power 
units, your value of Fq is also in absolute units. If you cannot measure P in absolute units, 
you can still measure a relative intensity profile in frequency and angle, normalized to the 
peak intensity. The resulting temperature is also in relative units.) 



VII. OTHER PHYSICAL CONSEQUENCES 

Nonequilibrium radiation temperature and its large magnitude in lasers have important 
theoretical and practical consequences. 

Consider the entropy production density of radiative transfer (matter-radiation coupling). 
The matter interacting with the radiation has a local temperature T. This entropy produc- 
tion density can be related to the radiation energy specific intensity Ii, and a entropy specific 
intensity Jin- 

- ""^J"^ + Sr + VB.= [ dudn^l - (k • V/, + -l)^ + (k ■ vj, + -X)} 

1 J c I c 

^/...0,{(k.v/„ + i4)(f -i)} 

= jd^d{l^{ (k . V4 + (-L- - i) } . (36) 

H and F are vector densities of radiative entropy and energy flux, respectively, and u and 
Sr are radiative energy and entropy densities, respectively. The radiation temperature T^{k) 
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is defined in Eqs. (|7j), ()13p . and ()19|). The first term on the first hne is the matter entropy 
production density; the second and third are the radiation entropy production density. 

Equation (j36|) is positive semidefinite and expresses the second law of thermodynamics 
in radiative transfer. The left-hand factor in the last integrand of Eq. (j36p is the energy 
transfer rate density, which is equal and opposite for matter and radiation.— For strictly 
thermal emission, this factor corresponds to the net cooling or heating of matter due to the 
radiation beam in the direction k. 

• If T > Ti,(k), it is positive: the matter cools, the field gains energy, and the entropy 
production is positive. 

• If the temperature inequality is reversed, the sign of this factor also reverses: the field 
loses energy, the matter is heated, and the entropy production remains positive. 

• The product in the integral is zero only when the matter and radiation temperatures 
are equal and no transfer occurs. 

The two other ways matter can transfer heat to other matter are by convection and diffu- 
sion. These mechanisms require matter-matter contact and gradients in intensive variables 
such as temperature. Radiative energy transport, on the other hand, depends only on the 
difference of the local matter and radiation temperatures at a single point in space. Even 
a low-power laser is very effective at transferring energy, because its radiation tempera- 
ture is so large. Energy transfer dominated by radiative processes is minimally affected by 
convection and conduction if the matter temperature is much smaller than the radiation 
temperature. 

Practical consequences of these results include the following. Because of the high beam 
temperature, radiative transfer dominates over the other, more destructive and undesirable 
heat transfer processes of convection and conduction, making lasers effective tools for surgery. 
There is neither the time nor the energy to induce gradients large enough to make convection 
or conduction important. The overall deposited energy is small, and the cut tissue suffers 
little damage. The cooling of atoms with lasers is also interesting in light of the results of this 
paper. In Sec. II V CI two cases of vanishing entropy were considered. Case 2 represents the 
laser radiation cooling the atoms, while Case 1 conforms to the behavior of the atoms being 
cooled. When the atoms decrease in entropy, according to Case 1, their temperature and 
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energy must also go down. But the radiation behaves according to Case 2: its temperature 
also goes down, but its entropy increases instead, while its power stays constant. 

The classical field limit of T^{^) is obtained by letting /i — in Eq. (jl3p and holding 
everything else constant. The result, A;BTj,(k) = hun = c^/,^(k)/z/^, is identical to the 
high-intensity or low- frequency limit, independent of h if written in terms of /^(k). The 
magnitude of the flux density |F| = F is 

F = jdudni,P{v-k)kBmk). (37) 

This result can be compared to the antenna or noise temperature familiar in radio-frequency 
electromagnetism.— In a one-dimensional system, the power Pnoise of a pure noise signal is 
associated with a temperature Tnoise: 

/"OO 

Pnoise = / dvkBT^o\se{v). (38) 

The result ()37|) is the three-dimensional, solid-geometry analogue of the radio-frequency 
antenna or circuit result (j38|) . It should be stressed that the antenna temperature is the 
temperature of the radiation, not the antenna material; just as Ti/(k) is not the temperature 
of the lasing gas, but of the laser light. 

The arguments of Sec. IVlmake clear that the noise temperature deflned in Eqs. fl37p and 
(|HHj) requires a randomized phase distribution for the noise signal, in keeping with the usual 
intuitive deflnition of noise. 

The results of this paper demonstrate some of the physical implications of temperature 
for nonequilibrium systems, in particular for radiation and laser light. More consequences 
can be inferred by applying the techniques and results presented here to other aspects of 
radiation, radiation transport, and lasers, and can be explored starting with the references. 
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FIG. 1: The gain function G{z) relating the radiation intensity, frequency, and temperature. 
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He-Ne Laser temperature 
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FIG. 2: The helium-neon laser radiation temperature T^iy, 9) for D = 1 mm, Az//? = 0.9 GHz, and 
Po = 1.0 mW. 
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He-Ne Laser Temperatures for Different Apertures 
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FIG. 3: The helium-neon laser radiation temperature T(z^Oi ^) for three aperture diameters D. The 
intensity/temperature profile narrows as D increases and the beam becomes more parallel. 
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